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1 . . PROJECT  ABSTRACT  . 

The  aim  of  the  proposed  research  is  the  performance  of  quantum 
mechanical  calculations  of  an  exact  and  approximate  nature  of  the  cross 
sections  and  rate  constants  of  elementary  chemical  reactions.  These  calcu¬ 
lations  are  of  importance  for  the  fundamental  understanding  of  the  nature  of 
chemical  reactions  and  for  advanced  technologies  of  interest  to  the  United 
States  Air  Force,  such  as  high-energy  chemical  lasers,  plume  technologv. 
and  the  chemical  effects  of  high-energy  radiation.  The  calculations  include 
processes  involving  three  atoms,  such  as  exchange  reactions,  collision-induced 
dissociation,  three-body  recombination,  electronic  branching  ratios,  tunneling, 
and  vibrational  quenching  by  reaction.  The  goal  of  this  work  is  to  develop  a 
detailed  understanding  of  such  processes  and  to  generate  benchmark  results 
which  can  serve  as  a  basis  for  testing  the  validity  of  approximate  methods. 
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2 .  SUMMARY . OF  PROGR E SS 

We  summarize  here  the  progress  of  the  last  year  in  the  work  supported 
by  this  contract. 

2. 1  Quantum  Mechanical  Calculations  of  Three-Dimensional  Reactive 

Scattering 

We  have  previously  developed*  the  conceptual  and  numerical  method¬ 
ology  for  performing  accurate  three-dimensional  quantum -mechanical  calcula¬ 
tions  of  the  corss  sections  of  triatomic  exchange  reactions  of  the  type  A  +  BC  ~ * 

2 

AB  +  C,  AC  +  B  and  applied  it  to  the  H  +  H2  system.  These  are  the  only 
accurate  three-dimensional  differential  reaction  cross  sections  published  so 
far.  A  modified  version  of  the  code  was  developed  for  systems  of  the  type  A2B, 
in  which  one  of  the  three  atoms  could  be  distinct  from  the  other  two.  During 
the  last  year  we  have  performed  extensive  test  calculations  on  the  H2D  system, 
i.  e.  ,  on  the  reaction  H  +  DH  H2  +  D,  HD  +  H,  and  studied  in  detail  the  num¬ 
ber  of  rotational  and  vibrational  basis  sets  and  the  matching  program  param¬ 
eters  required  for  achieving  convergence  of  the  results.  This  testing  work 
was  done  on  our  VAX  11/780  computer.  The  production  calculations  that  must 
be  performed  in  the  future  are  two  lengthy  for  such  a  computer,  and,  as  a 
result,  we  are  transferring  them  to  an  IBM  3032  machine  available  for  this 
purpose.  This  transfer  required  further  modifications  of  the  code,  as  the 
size  of  its  fast  memory  is  smaller  than  that  of  the  VAX.  The  necessary  modifi¬ 
cations  have  been  made  and  are  currently  being  tested.  We  expect  to  reach 
production  stage  on  the  IBM  computer  for  zero  total  angular  momentum  in  the 
next  few  months.  The  process  is  lengthy  and  painstaking  in  view  of  the  size 
and  complexity  of  the  codes. 
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2.  2  The  Use  of  Hyperspherical  Coordinates  for  Collinear  Exchange 

Reactions 

We  had  previously  developed  a  new  method  for  performing  quantum 

mechanical  calculations  on  collinear  triatomic  collisions.  Over  the  last  year 

this  method  has  been  extensively  tested  on  the  H  +  H>  and  F  +  H2  systems,  and 

found  to  be  significantly  more  efficient  than  previous  ones.  A  description  of  it 

3 

and  new  results  on  the  H  +  system  have  been  published,  and  a  reprint  is 
appended.  In  particular,  it  requires  few^r  vibrational  channels  for  converg.-v.c0, 
compared  with  other  methods.  For  the  two  systems  just  mentioned,  about  1/2 
the  usual  number  is  sufficient,  and  for  a  heavy-light- heavy  type  system,  such 
as  I  +  HI,  1/6  of  the  usual  number  does  the  job.  Since  the  time  required  for 
such  calculations  typically  increases  with  the  cube  of  the  number  of  coupled 
channels,  the  new  method  achieves  time  economies  of  23  or  more,  which  is  a 
substantial  improvement. 

For  heavy -light -heavy  systems,  such  as  I  +  HI,  I  +  HC1,  etc. ,  the  usual 
methods  not  only  require  large  numbers  of  coupled  channels  (40  or  more),  but 
also  have  great  difficulty  in  yielding  reliable  results.  The  central  reason  is  the 
smallness  of  the  skew  angle  of  the  corresponding  potential  energy  surfaces, 
which  is  only  7°  for  the  I  +  HI  system  (as  compared  with  60°  for  H  +  Kj).  The 
hyper  spherical  coordinate  method  has  no  special  difficulties  with  such  a  system 
and  has  yielded  high  quality  results  for  I  +  HI  using  only  six  coupled  channels. 

4 

The  corresponding  results  and  discussion  have  been  published  recently,  and 
a  reprint  is  appended. 

The  method  of  hyperspherical  coordinates  permits  calculations  to  be 
done  at  high  energies  without  deterioration  of  the  quality  of  the  results.  We 
have  calculated  from  such  results  deactivation  rates  for  vibrationally  excited 
Hj  (v)  molecules  by  reactive  and  nonreactive  collinear  collisions  with  H  atoms 
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over  the  temperature  range  200°  K  to  2000"  K  for  v  =  1  through  5.  The  corres¬ 
ponding  Arrhenius  plots  show  noticeable  curvature  at  high  temperatures  and 
obey  an  exponential  gap  law  quite  well. 

5  6 

A  system  of  considerable  experimental3  and  theoretical0  interest  is  the 
Cl  +  HC1  ~~ 1 "  C1H  +  Cl,  another  heavy -light- heavy  system,  with  a  skew  angle  of 
13.  5°.  We  have  performed  reaction  probability  and  rate  constant  calculations 
on  this  system  on  two  LEPS  surfaces,  with  barrier  heights  of  6.  2  and  1.  9  kcal/ 
mole,  respectively.  We  found  that  the  rate  constant  for  reaction  of  HC1  (v  =  1) 
is  higher  than  that  of  HC1  (v  =  0)  for  the  larger  barrier  height  surface,  but  that 
the  reverse  is  true  for  the  lower  one.  In  fact,  the  HC1  (v  -  1)  rate  on  the  low 
barrier  surface  is  smaller  than  on  the  high  barrier  one.  This  is  an  unusual 
result,  related  to  the  heavy-light-heavy  nature  of  the  system. 

Reactions  of  the  general  type  X  +  HY(v)  -*  XH  +  Y,  where  X  and  Y  are 

7 

halogens,  are  of  substantial  experimental  interest.  We  have  also  performed 
some  preliminary  calculations  on  the  Cl  +  HBr  system,  which  has  a  skew  angle 
of  11.  5°.  Convergence  to  about  1%  in  the  translational  energy  range  0  to  0.  4  eV 
was  achieved  with  six  basis  functions  only,  whereas  the  use  of  a  previous 

8 

method  based  on  the  one  developed  in  our  laboratory  required  40  basis  functions, 
an  additional  demonstration  of  the  superior  efficiency  of  the  hyperspherical 
method. 

2.  3  The  Use  of  Hyperspherical  Coordinates  for  Collinear  Collision- Induced 

Dissociation 

We  had  previously  developed  a  method  for  rigorously  discretizing  the 
dissociation  continuum  in  collision-induced  dissociation  and  three-body  recom¬ 
bination  processes.  During  the  last  year  we  obtained  accurately  converged 
calculations  for  the  collision- induced  dissociation  probabilities  on  a  model 


symmetric  collinear  potential  energy  surface.  The  method  has  functioned 
remarkably  well,  the  probabilities  having  converged  to  within  ±  0.  02  for 
reagent  translation  energies  up  to  three  times  the  dissociation  energy  with 


only  eight  channels  of  each  parity  used  to  represent  the  continuum.  The  results 

9 

are  in  press,  and  a  preprint  is  appended. 

2.  4  Collinear  Reaction  Probabilities  for  the  H  +  HI  ~ *  E,  +  I  Reaction 

This  reaction  is  of  considerable  interest  because  its  exothermicity  is 
sufficient  to  produce  electronically  excited  I  atoms.  As  a  preliminary  step  in 
its  study,  we  have  done  calculations  on  a  single  electronically  adiabatic  poten¬ 
tial  energy  surface,  preparatory  to  a  two-surface  investigation.  In  a  previous 
study  of  this  system,  a  potential  energy  surface  having  unrealistic  wells  in 
the  reagent  and  product  channel  regions  was  used  and  resulted  in  pronounced 
sharp  resonant  structure.  To  avoid  this,  we  constructed  a  more  realistic 
surface,  of  the  rots'-’ng  Morse-cubic  spline  type,  **  devoid  of  such  unphysical 
wells.  The  sharp  resonances  disappeared,  as  expected,  and  it  was  found  that 
the  bulk  of  the  ^  product  was  formed  in  the  v  =  2  state.  This  finding  may  be 
a  useful  guide  for  the  deisgn  of  crossed  beam  experiments  on  this  important 
system. 

2.5  Three-Dimensional  Reactive  Scattering  Calculations  Using  Symmetrized 

Hyperspherical  Coordinates 

Over  the  last  year  substantial  progress  was  achieved  in  the  application 
of  the  hyperspherical  coordinate  method  to  the  three-dimensional  H  +  H, 
reaction.  The  changes  indicated  in  the  last  progress  report,  involving  a  modifi¬ 
cation  of  the  method  used  for  calculating  the  surface  functions,  were  imple¬ 
mented.  They  resulted  in  conservation  of  flux  for  zero  angular  momentum 
quantum  number  of  1%  or  better,  converged  with  respect  to  the  number  of 
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basis  functions.  We  find  that  ten  such  channels  are  sufficient  over  the  C.  3  to 

o 

0.  7  eV  total  energy  range,  whereas  previously  30  to  40  channels  were  necessary. 
This  is  a  significant  improvement  of  efficiency,  which  can  result  in  a  decrease 
of  computer  time  by  a  factor  in  excess  of  10.  However,  the  results  differ  by 
about  10%  from  those  obtained  previously.  We  believe  that  this  difference  is 
due  to  a  deterioration  of  the  quality  of  the  surface  functions  for  hyperdistances 
in  excess  of  5  bohr,  and  we  are  currently  implementing  a  change  in  this  calcu¬ 
lation  for  such  distances  aimed  at  correcting  this  defect.  It  is  expected  that 
the  code  will  become  operational  in  the  not  too  distant  future. 
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DESCRIPTION  OF  PROPOSED  RESEARCH 


We  intent  to  continue  in  the  next  year  or  two  to  do  research  along  the 

Lines  described  in  the  proposal  dated  31  August  1976,  as  follows: 

1.  Continue  the  three-dimensional  calculation  on  the  D  +  YL,  system  in 
order  to  generate  differential  and  integral  cross  sections  and  rate 
constants,  comparing  them  with  those  for  H  +  Ho,  to  obtain  accurate 
isotope  effects  on  the  same  potential  energy  surface.  Also,  we  will 
test  the  validity  of  isotope  effect  theories. 

2.  Perform  three-body  recombination  calculations  on  the  same  collinear 
surface  for  which  collision-induced  dissociation  calculations  were 
performed  last  year.  Also,  study  the  partitioning  of  reagent  transla¬ 
tional  energy  among  the  products  of  collision-induced  dissociation. 

3.  Do  a  systematic  study  on  the  dynamic  characteristics  of  the  X  +  HY  — 
XH  -r  Y  collinear  reactions,  where  X  and  Y  arc  halogen  atoms,  as  a 
function  of  features  of  the  potential  energy  surface,  using  the  new 
hyperspherical  coordinate  method. 

4.  Study  the  electronic  branching  ratio  for  the  collinear  H  +  HI  ~ *  H,  +  I 
(2I>3/2,2F>i  /2  )  reaetion  for  m°del  coupling  surfaces,  to  elucidate  the 
conditions  under  which  the  production  of  electronically  excited  products 
is  formed. 

5.  Continue  the  development  of  the  three-dimensional  hyperspherical 
coordinate  code  and  apply  it  to  the  H  +  H2  and  D  +  Hj  reactions  as 
benchmark  systems.  Also,  extend  it  to  the  T  +  HD  one.  Evaluate  its 


convergence  with  respect  to  size  of  basis  set  and  its  computational 
efficiency. 

6.  Perform  ^-conserving  three-dimensional  approximate  calculations 

1  2 

based  on  the  older  arrangement  channel  method  ’  and  the  newer 


h^erspherical  method,  as  well  as  corresponding  first-order  perturba¬ 
tion  calculations  whose  reference  potential  contains  the  important 
dynamic  features  of  the  correct  potential. 

7.  Perform  accurate  three-dimensional  calculations  on  the  Cl  +  ¥1,  reaction 
at  low  energies  by  either  the  arrangement  channel  or  hyperspherical 
coordinate  methods. 

8.  Perform  calculations  on  streamlines  of  probability  current  density  and 
tunneling  fractions  in  order  to  assess  the  reaction  pathways  on  the 
potential  energy  surface  for  the  reactions  mentioned  in  3  above. 

The  order  in  which  these  several  projects  will  be  pursued  will  depend 
on  available  manpower  and  relative  progress,  but  it  is  expected  that  a  signifi¬ 
cant  effort  will  be  carried  on  on  most,  if  not  all,  of  them. 
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4^ _ 

The  Project  Director  will  continue  to  be  Professor  Aron  Kuppermann, 
whose  curriculum  vitae  and  publications  list  were  included  in  the  proposal 
dated  31  August  1976.  The  following  additional  personnel  are  currently  work¬ 
ing  on  different  facets  of  the  project. 

Dr.  Joseph  K.  C.  Wong,  postdoctoral  fellow,  full-time 

Dr.  Nancy  M.  Harvey,  postdoctoral  fellow,  1/3  time 

Mr.  Jack  A.  Kaye,  graduate  student 

Ms.  Diane  M.  Hood,  graduate  student 

Mr.  Huey-Jenn  Chiang,  graduate  student 

Mr.  Loudon  L,  Campbell,  undergraduate  student. 
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5.  BUDGET  (A FOSR):  1  July  1C81  through  30  June  1982 
Salaries 


Aron  Kuppermann,  Principal  Investigator 

One  Postdoctoral  Research  Fellow,  full-time, 

12  months  @  $13,  500/year 

Two  Graduate  Research  Assistants,  full-time 
12  months  @$ 6060/year 

Secretary,  12  months,  1/4  time 

Total  Salaries 

Staff  Benefits 

24.  6%  of  Total  Salaries 
Supplies  and  Expense 

Stationery  supplies,  photocopying,  graphic  arts,  illustrations 
Page  charges  and  reprint  costs  for  publications 

Computing  Costs 

700  hours  @  $  14/hour  (DEC  VAX  11/780) 

Travel  -  Domestic 

Meetings  and  Conferences 

Total  Direct  Costs 

Overhead 

48%  of  Total  Direct  Costs 

TOTAL  COSTS 


$  5,600 
13,500 

12,120 

4.800 

$36,020 

8,861 

1,700 

1,200 

9.800 

1,500 

$59,081 

28, 359 
$87,440 
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I-.-.  LIST  OF  PUBLICATIONS 

The  following  are  publications  describing  work  performed  under  this 
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